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ABSTRACT 

Riboswitches are a newly discovered large family of 
structured functional RNA elements that specifically 
bind small molecule targets out of a myriad of 
cellular metabolites to modulate gene expression. 
Structural studies of ligand-bound riboswitches by 
X-ray crystallography and nuclear magnetic reson- 
ance (NMR) spectroscopy have provided insights 
into detailed RNA-ligand recognition and inter- 
actions. However, the structures of ligand-free 
riboswitches remain poorly characterized. In this 
study, we have used a variety of biochemical, bio- 
physical and computational techniques including 
small-angle X-ray scattering and NMR spectroscopy 
to characterize the ligand-free and ligand-bound 
forms of SAM-II riboswitch. Our data demonstrate 
that the RNA adopts multiple conformations along 
its folding pathway and suggest that the RNA 
undergoes marked conformational changes upon 
Mg 2+ compaction and S-adenosylmethionine (SAM) 
metabolite binding. Further studies indicated that 
Mg 2+ ion is not essential for the ligand binding but 
can stabilize the complex by facilitating loop/stem 
interactions. In the presence of millimolar concen- 
tration of Mg 2+ ion, the RNA samples a more 
compact conformation. This conformation is near 
to, but distinct from, the native fold and competent 
to bind the metabolite. We conclude that the forma- 
tion of various secondary and tertiary structural 
elements, including a pseudoknot, occur to seques- 
ter the putative Shine-Dalgarno sequence of the 
RNA only after metabolite binding. 

INTRODUCTION 

S-Adenosylmethionine (SAM), as a major methyl donor, 
plays a central role in signal transduction and transcrip- 
tional regulation of many critical processes in sulfur 



metabolism of certain bacteria (1). While such regulatory 
mechanisms of feedback inhibition performed by proteins 
are sophisticated and have long been recognized (2,3), it is 
becoming increasingly clear that non-protein coding 
RNAs, such as riboswitches, also play an important role 
in the regulation of gene expression (4-8). 

Riboswitches are structured functional RNA elements 
that are typically located at the 5' non-coding regions of 
messenger RNAs (9) and comprise an aptamer domain 
and an expression platform domain. The aptamer 
domain functions as a sensor to recognize and specifically 
bind small molecule targets out of a myriad of intracellu- 
lar metabolites, whereas the expression platform modu- 
lates gene expression in response to the binding events in 
the aptamer domain. It is thought that a change in con- 
formation can lead to either inhibition of transcription 
elongation via formation of terminator structure or 
blockage of translation through sequestering of the 
ribosome binding site (4,10,11). This presumed regulatory 
strategy is likely achieved through the interplay of two 
mutually exclusive secondary structures of the expression 
platforms that lead to either expression or repression of 
the downstream genes encoding for the enzymes respon- 
sible for the transport or biosynthesis of the target metab- 
olites (12,13). 

Although significant progress has been made in under- 
standing how riboswitches recognize their ligands (14-23), 
the nature of the unbound state is still not very well 
defined (24-28). Instead, several conflicting perspectives 
about the unbound state have arisen (27). First, both the 
apo and the bound states of a Thermotoga maritima asd 
lysine riboswitch have identical structures with the ligand 
binding pocket completely buried within the molecule, 
making it hard to rationalize how the ligand interacts 
with the riboswitch and enters into the binding pocket 
(22,29,30). In addition, nuclear magnetic resonance 
(NMR) spectroscopy and fluorescence studies of the 
Bacillus subtilis queC and Fusobacterium nucleatum 
preQl riboswitches have revealed that the receptors are 
unstructured in the absence of ligand binding and that 
addition of preQl induced secondary structural 
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rearrangement and global changes in the architecture 
of the riboswitches to shut down gene transcription 
(21,31). Ligand-induced folding has been investigated 
thermodynamically and kinetically, but it is still unclear 
how the ligand can recognize such a disorganized bind- 
ing pocket (31). Moreover, for translational acting 
riboswitches such as S M k box (SAM-III), distinct 
mutually exclusive RNA conformations are differentially 
populated in the presence or absence of the metabolite. It 
remains unclear how the metabolite triggers a switch from 
a different structure to a complex and highly specific 
binding pocket (32-34). 

To better understand the functional mechanism of a 
riboswitch and its folding pathway in the unliganded 
state, we analyzed a SAM-II riboswitch with a sequence 
located upstream of the metX gene (encoding homoserine 
acetyltransferase) from the Sargasso Sea metagenome 
(35). X-ray crystal structural analysis indicated that 
three helical elements, two loops and a classic H-type 
pseudoknot all come together to form an intricate 
tertiary structure with the SAM binding pocket located 
in the center of helix P2b (Figure 1) (23). Striking 
features of the SAM-bound RNA structure are that the 
helical segments stack upon each other to create an overall 
extended conformation and that the 3'-end corresponding 
to the Shine-Dalgarno (SD) sequence is stabilized by 
base-pairing and tertiary interactions with loop 1 (23). 
These, together with the RNAs binding selectivity and 
the fact that SAM-II is one of the only riboswitches con- 
taining the sequences of both the aptamer and the expres- 
sion platform, make SAM-II an attractive system to 




Figure 1. Secondary structure and X-ray crystal structure of SAM-II 
RNA (Protein Data Bank ID: 2qwy). (a) Metabolite-bound secondary 
structure of SAM-II riboswitch proposed from the X-ray crystal struc- 
ture. The RNA used in this study has two single point mutations (U1G 
for improving transcription efficiency and G52A for enhancing 3'-end 
cleavage by the 3' H8V ribozyme) compared with its wild-type. The 
X-ray structure contains two additional double point mutations (C6G 
U7C and A25G G26U) for adding a cesium/iridium ion binding site for 
phasing (23). These are indicated by asterisk, (b) X-ray crystal structure 
of SAM-II riboswitch with the same color coding as the secondary 
structure (Protein Data Bank ID: 2qwy) (23). 



investigate the dynamics and alternative conformational 
switching upon ion and SAM binding. 

We hypothesized that the predicted ribosome binding 
site of the free form RNA is exposed and accessible to 
the ribosome for recognition and that the ligand-free 
SAM-II RNA adopts multiple conformations that are 
modulated by ions and ligands. In this study, we charac- 
terize the SAM-II riboswitch under different folding con- 
ditions using small-angle X-ray scattering (SAXS) 
combined with NMR spectroscopy. These data suggest 
that the RNA undergoes marked conformational 
changes upon Mg 2+ compaction and SAM metabolite 
binding and that Mg 2+ binding induces a state that is 
competent to bind metabolite. Only after metabolite 
binding does the formation of various secondary and 
tertiary structural elements including a pseudoknot 
occur to sequester the SD sequence of the RNA. In the 
course of this work, a similar study using fluorescence and 
NMR spectroscopy and single molecule fluorescence 
energy transfer (smFRET) imaging revealed a dynamic 
unliganded SAM-II riboswitch (36). Our work comple- 
ments and extends this recent study in two important 
aspects: SAXS provides direct evidence for at least four 
conformational states of varying compactions; NMR 
analysis provides greater insights into the nature of the 
secondary and tertiary structural rearrangement of the 
RNA in the presence of Mg 2+ ions and metabolite. 

MATERIALS AND METHODS 

Reagents and materials 

Oligonucleotides for the polymerase chain reaction were 
purchased from Integrated DNA Technologies, Inc. The 
SAM was obtained from Sigma and Tag DNA polymer- 
ase from New England Biolabs, Inc. A plasmid containing 
the SAM-II RNA sequence along with a self-cleaving 
transacting H8V ribozyme at the 3'-end was a gift from 
Dr Robert T. Batey of the University of Colorado at 
Boulder. 

RNA preparation 

The RNA used in this study is a variation of the wild type 
RNA described by Gilbert et al. (23), in which two single 
point mutations were introduced to change uracil 1 to 
guanine for improving transcription efficiency and 
guanine 52 to adenine for enhancing 3'-end cleavage by 
the 3' H8V ribozyme. The RNA was transcribed as 
described. 

Size-exclusion chromatography 

Chromatographic separation was performed with the 
AKTA™ Design system with a detector at a wavelength 
of 254 nm. The conformational homogeneity of the RNA 
sample was checked using size exclusion chromatography 
(SEC) with a Superdex 200 HR 10/30 column (GE 
Healthcare; column height: 30 cm, column diameter: 
1.0 cm, column volume: 23.562 ml). Samples of SAM-II 
riboswitch in 50 mM sodium— 2 -(A^-morpholino)etha- 
nesulfonic acid (Na-MES), pH 6.2, 50 mM NaCl were 
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denatured at 90° C for 3min, snap cooled and incubated 
on ice for lOmin, then Mg 2+ and/or the metabolite of 
varying concentrations were added, and the mixtures 
were incubated on ice for 30min before loading onto the 
column. A flow rate of 0.5ml/min was used in all the 
chromatographic steps performed at 4°C. 

Electrophoretic mobility shift assay 

Samples of SAM-II riboswitch used in electrophoretic 
mobility shift assay were prepared as described above. 
Electrophoresis was performed using 12% native poly- 
acrylamide gel electrophoresis (PAGE) in 90 mM Tris, 
90 mM MES, pH 8.0, in the absence or presence of 
50 mM NaCl and varying concentrations of MgCl 2 or 
the metabolite at 4°C at 300 V. The RNA samples and 
gel matrix contained identical concentrations of MgCl 2 
and the metabolite, and running buffer was kept 
recirculating every 1 h. The gels were stained with 
ethidium bromide and images were scanned on a 
STORM phospholmager (GE Healthcare). 

Small angle X-ray scattering experiments and 
data analysis 

Both SAXS and wide angle X-ray scattering (WAXS) data 
were collected at beam line 12-ID of the Advanced Photon 
Source at the Argonne National Laboratory. The wave- 
length, A, of the X-ray radiation was set to 1.033 A. 
Detailed protocols of SAXS measurement have been 
described elsewhere (37). In brief, the SAM-II riboswitch 
RNA was dissolved in buffer A (50mM Na-MES at 
pH 6.2). Samples were heated at 90°C for 3min followed 
by cooling on ice for lOmin. MgCl 2 and/or SAM ligand 
were added to prepare four RNA sample buffer condi- 
tions: Al (buffer A without MgCl 2 or ligand), A2 
(buffer A with 10 mM MgCT), A3 (buffer A with 
0.5 mM ligand) and A4 (buffer A with lOmM MgCl 2 
and 0.5 mM ligand). To ensure identical buffer conditions, 
each sample was then exchanged with the appropriate 
buffers (A1-A4) six times using an Amicon 
microconcentrator filter (Millipore). RNA samples were 
prepared in three concentrations (0.8, 1.8 and 3.6mg/ml) 
for each of the four buffer conditions before shipping on 
dry ice, therefore, making a total of 12 samples. 

The g-range used for data analysis was between 0.006 
and 0.300 A -1 . The linearity of the Guinier plot was used 
to evaluate the quality of the SAXS data. WAXS data 
were used to guide accurate background subtraction for 
the corresponding SAXS data (38). A total of 20 images 
were taken for each sample condition to get good statistics 
and converted to ID scattering profiles. The resulting scat- 
tering data sets were averaged before buffer background 
subtraction using IGOR PRO (WaveMetrics) software. 
The reported radius of gyration (R s ) values were derived 
from a Guinier approximation in the (/-range such that 
<7max x -R g < 1-3. The GNOM software (39) was used to 
determine R g values, which were in good agreement with 
those obtained from Guinier analyses. This software was 
also used to perform pair distance distribution function 
(PDDF) plots. 



Low resolution ah initio model reconstructions of 
SAM-II riboswitch 

GNOM (39), DAMMIN (40) and DAMAVER (41) from 
the ATSAS 2.4 suite of programs were used to perform 
ab initio model reconstructions of the RNA under the four 
solution conditions described above. To avoid underesti- 
mation of the molecular dimension, we chose the param- 
eter ^? max , the upper end of distance r, such that the 
resulting PDDF has a short, near-zero-value tail at large 
r. The maximum distances (D max ) were estimated to be 
112, 72, 70 and 68 A from the PDDF for the SAM-II 
RNA under the four experimental conditions (see 
'Results' section). We used the program DAMMIN to 
obtain an approximate molecular envelope. To avoid dis- 
tortion caused by possible underestimation of D max , we 
usually set the R ma x to be about 10 A greater than Z) max . 
All reconstructions for SAM-II under the four conditions 
were performed in the 'slow' mode. The resulting struc- 
tural models were averaged using the program package 
DAMAVER (41), which generates an average and 
filtered dummy atom model. In this program, the 
normalized spatial discrepancy (NSD) values between 
each pair of models were computed. The model with the 
lowest average NSD with respect to the rest of the models 
was chosen as the reference model. The remaining models 
were superimposed onto the reference model using 
SUPCOMB. Possible outliers identified by NSD criteria 
were discarded. A total of 20 independent DAMMIN runs 
were performed, and the resulting bead models were 
averaged using DAMAVER. CRYSOL (42) was used to 
calculate the R g value and the SAXS profile for the 
SAM-II crystal structure. 

NMR sample preparation and NMR experiments of 
SAM-II RNA 

The purified RNA was lyophilized and dissolved in water, 
and the corresponding amount of buffer stock solutions 
was added. The RNA sample solutions were refolded 
using identical conditions as described above. Then 
8-10% D 2 0 and various concentrations of SAM ligand 
and Mg 2+ were added and the sample volume was 
adjusted to 260 ul with ddH 2 0 in Shigemi NMR tubes. 
The final NMR buffer contains 10 mM sodium phosphate, 
pH 6.2 and 100 mM NaCl. The sample concentration was 
~0.5— 1 mM for NMR experiments as judged by the mag- 
nitude of UV absorbance at 260 nm and the calculated 
extinction coefficient (668 300 M -1 cm -1 ) of the individual 
RNAs base composition. 

NMR experiments were performed on a Bruker Avance 
600 MHz spectrometer equipped with actively shielded 
z-axis gradient triple resonance probe. NMR spectra 
were recorded at a temperature of 298 K using the 
WATERGATE water suppression pulses (43). 'H- 15 N 
HSQC, 'H- 13 C HSQC, 2 Jhn-'H-^ 5 N HSQC (44) and 
2D NOESY experiments were carried out using standard 
pulse sequences (45). The 2 J NN — HNN— COSY experi- 
ments were performed as described (46). Homonuclear 
NOESY spectra used to determine imino proton 
connectivities were recorded with a mixing time of 
300 ms, with 256 complex points in the indirect dimension. 
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'H- 15 N HSQC spectra were recorded with 256 complex 
points in the indirect dimension. All spectra were pro- 
cessed and analyzed using Bruker's TOPSPIN 2.1, 
NMRPipe/NMRDraw (47) and N MR View (48). 

RESULTS 

We first performed conformational heterogeneity analysis 
of the SAM-II RNA sample under different solution con- 
ditions using analytical SEC and electrophoretic mobility 
shift assay in the absence or presence of Mg 2+ and/or 
SAM metabolite. SAXS data were then collected on 
SAM-II riboswitch in the absence or presence of Mg 2+ 
and/or SAM metabolite. Based on the SAXS data, we 
reconstructed four ab initio low-resolution electron 
density models corresponding to four different folding 
states of the RNA. Finally, we characterized the free 
and native states of the RNA using NMR spectroscopy 
to complement the SAXS data. Models for the RNAs 
folding and regulatory mechanism were proposed based 
on our data. 

Conformational heterogeneity analysis of the RNA sample 

Preparation of homogenous macromolecule samples is 
critical for biochemical and biophysical studies. 
Generally, RNA synthesis involves several steps including 
an in vitro transcription with T7 RNA polymerase, puri- 
fication by a denaturing PAGE, dialysis and buffer 
exchange and storage at — 20°C, denaturation and 
renaturation before biochemical and biophysical studies. 
This production protocol, although robust and well- 
established, suffers from several disadvantages in 
addition to probable aggregation of RNA samples (49). 
First, RNA molecules derived from denaturing PAGE 
purification contain acrylamide oligomer contaminations 
that will interfere with the SAXS scattering analysis and 
compromise interpretation of NMR data (50,51). Second, 
RNA molecules, as a result of the renaturation process, 
may misfold into kinetically trapped species that can result 
in conformational heterogeneities (51-53). However, 
accurate SAXS data interpretation and model reconstruc- 
tion require homogenous samples with defined conformers 
in solution because the intensity of the scattering data 
results from population weighted averages of parameters 
from individual components (54). As a result, conform- 
ational heterogeneity of the RNA increases the possibility 
of misinterpretation of SAXS data, leading to reconstruc- 
tion of artif actual bead models. 

To avoid this problem, we examined the conformational 
homogeneity of the SAM-II RNA before SAXS analyses 
using SEC on fast performance liquid chromatography 
(FPLC). The RNA was first synthesized by in vitro tran- 
scription with T7 RNA polymerase and purified by 
denaturing PAGE to homogeneity. A folding protocol 
of snap-cooling on ice after denaturing at 90° C for 3min 
was used to test its effect on the RNA folding in the 
absence or presence of Mg 2+ ion and/or SAM metabolite. 
The elution profile of the SEC run shows a single peak of 
the SAM-II RNA in the absence of Mg 2+ and metabolite, 
indicating that no minor conformer exists or were too 



insignificant to be detected under our experimental condi- 
tions (Figure 2a). 

More analytical SEC analyses were conducted for the 
RNA in the presence of different concentrations of Mg 2+ 
ion and/or saturating amounts of metabolite to examine 
the conformational homogeneity of the RNA under these 
solution conditions. An overlay of SEC chromatograms of 
the RNA in the presence of different Mg 2+ ion concentra- 
tions as a function of retention time reveals that no 
shoulder peak is present, an indication that the RNA is 
a monomer and homogenous under all of the solution 
conditions (Figure 2b). In the absence of Mg 2+ , the 
SAM-II riboswitch elutes from the column at ~14.5ml. 
Addition of 1 mM Mg 2+ ion to the mobile phase resulted 
in a slight delay in the retention time, suggesting a 
decrease of particle hydrodynamic radius. As the Mg 2+ 
concentration increases, the RNA elutes with longer re- 
tention time from the column. Compaction of SAM-II 
riboswitch as a function of Mg 2+ concentration reveals a 
folding transition with Mg 1/2 ~6mM (Figure 2c). Similar 
compaction was also observed for the RNA in the 
presence of saturating amount of the metabolite 
(Figure 2d), monovalent ion or other divalent ion such 
as Ca 2+ (Supplementary Figure SI). This is most likely 
because the Mg 2+ or metabolite induced compaction of 
the RNA caused a decrease in particle size, making the 
RNA travel a more circuitous path through the SEC sta- 
tionary phase, and as a result exhibited a delay in reten- 
tion time. Given the different behaviors of RNA and 
protein on SEC, this delay suggested remarkable 
changes of the particle size or hydrodynamic radii, 
although there is no standard curve for single-stranded 
RNA on a Superdex 200 (10/30) column to quantitatively 
monitor the change of the hydrodynamic radius (55). As a 
second method, the RNA was analyzed on native PAGE 
under a series of Mg 2+ ion and/or the metabolite concen- 
trations to examine their effects on the folding and con- 
formational change of the RNA. 

Native gel mobility shift assay 

Native PAGE is a simple but useful technique to test RNA 
folding and conformational changes at a gross level, and 
therefore applicable to riboswitches that change conform- 
ation upon binding metal ions or specific metabolites. If 
an alternative state of the RNA is present after refolding 
and the exchange between the two states is slow relative to 
the migration rate of the RNA through the gel, then bands 
with different migration rates will be observed, corres- 
ponding to different conformational states; otherwise, a 
single band of average mobility will be observed (56). As 
reported by Batey's group, the binding affinity between 
the SAM-II riboswitch and the metabolite is high (in 
the low micromolar range) in the presence of magnesium 
ions, indicating that once bound to the metabolite, the 
RNA likely adopts a compact conformation that should 
migrate more rapidly through the gel than the unbound 
RNA (57). 

An estimate of the relative migration rate of the SAM-II 
RNA was obtained by comparing the native gel electro- 
phoretic mobilities in the absence or presence of Mg 2+ or 
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Figure 2. Analysis of the homogeneity and conformation of the SAM-II RNA. (a) SEC elution profile for SAM-II riboswitch in the absence of Mg 
ion and the metabolite, (b) SEC elution profiles for SAM-II riboswitch in the presence of different concentrations of Mg 2+ ion (mM) in the mobile 
phase. Increasing MgCl 2 concentration resulted in preferential retention of the SAM-II riboswitch on the column, presumably because of decreasing 
hydrodynamic radius, (c) Compaction of SAM-II riboswitch as a function of MgCl 2 concentration reveals a folding transition with Mgi /2 = 6mM. 
(d) SEC elution profiles for SAM-II riboswitch in the presence of different molar ratios of the metabolite to SAM-II riboswitch. Addition of the 
metabolite caused a delay in the retention time of the RNA, consistent with a decrease in hydrodynamic radius of the RNA. Vertical axes for (a), (b) 
and (d) are in mAU at 254 nm. 



different stoichiometric amounts of the SAM metabolite 
(Figure 3). In the absence of Mg 2+ ion and metabolite, the 
RNA exhibits a homogenous single band. Addition of dif- 
ferent stoichiometric amounts of the metabolite, in the 
absence of Mg 2+ ion, did not change the electrophoretic 
mobility of the RNA (Figure 3a). In the presence of 
0.2 mM Mg 2+ without the metabolite, the RNA migrates 
as a homogenous band again. However, in the presence of 
unsaturating metabolite, the RNA adopts two conform- 
ations: an unbound form and a compact conformation. 
Upon adding saturating amount of the metabolite, the 
RNA switches completely to the compact conformation 
and the relative electrophoretic mobility of the compact 
complex is substantially accelerated. Consequently, it 
migrates more rapidly compared with the unbound 
RNA molecules. In the presence of 1 mM Mg 2+ and dif- 
ferent stoichiometric amounts of the metabolite, the RNA 
exhibits similar migration patterns to that observed in the 



presence of 0.2 mM Mg (Figure 3a), suggesting that the 
metabolite plays a critical role in RNA conformational 
change in the presence of low Mg 2+ ion concentration. 
Using bromophenol blue as a standard, we compared elec- 
trophoretic mobilities of different bands to that of the dye 
to obtain relative mobility (R { ) values. Notably, the Rf 
values of the RNA are different in the absence or 
presence of different Mg 2+ ion concentrations, indicating 
that Mg 2+ ion indeed plays a key role in promoting com- 
paction of the RNA molecule. Surprisingly, in the 
presence of 10 mM Mg 2+ , the relative electrophoretic 
mobilities of the RNA bound to different stoichiometric 
amounts of metabolite are nearly identical (Figure 3a). 
These observations suggest that at high Mg 2+ concentra- 
tions, the RNA molecule is already compacted to such an 
extent that addition of the metabolite has only an incre- 
mental effect on further folding the RNA, consistent with 
the results from SEC analyses. 
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Figure 3. Native PAGE of SAM-II riboswitch in 90mM Tris-MES, pH 8.0 in the absence (a) and presence (b) of 50mM NaCl. 2/1, 1/1, 1/5 
and 1/10 indicates the molar ratios between the RNA and the SAM metabolite. 



The mobility shift pattern is slightly different in the 
presence of 50 mM Na + ion. The RNA shifts negligibly 
upon addition of the metabolite in the absence of Mg~ + 
ion. In the presence of 1 mM Mg 2+ ion, the 
metabolite-induced migration shift is not as dramatic as 
that observed in the absence of 50 mM monovalent 
(Figure 3b). Taken together, these data suggest that sub- 
stantial rearrangement of the RNA occurs upon binding 
saturating amount of the metabolite or high concentration 
Mg 2+ ion and reveal that the RNA adopts different con- 
formations in the absence or presence of Mg 2+ and 
saturating amount of metabolite. The native gel shift 
assay provides a hint about big conformational changes. 
Further studies to address and provide models of the 
riboswitch folding as a function of Mg 2+ and the metab- 
olite by using SAXS and NMR analysis, however, indicate 
that there are still some conformational differences of the 
RNA under the conditions. 

Following the conformational changes of the SAM-II 
riboswitch induced by the metabolite and Mg 2+ ion 
with SAXS 

Solution scattering profiles were collected on SAM-II 
riboswitch in the absence or presence of 10 mM Mg 2+ 
and 0.5 mM SAM metabolite after passing the samples 
through a gel filtration column to ensure that the RNAs 
were monomeric and conformationally homogeneous. 
Four scattering curves representing the four buffer condi- 
tions described earlier (A1-A4) that were overlaid 
exhibited significant differences, indicating the existence 
of four possible conformational states along the 
RNA-folding pathway (Figure 4a). The SAXS data were 
analyzed using the linearity of the Guinier plot (a plot of 



natural log intensity [log 1(g)] versus q~, where q is the mo- 
mentum transfer). From the scattering data, the R g was 
calculated at low q values in the range of qR g < 1.3, using 
the Guinier approximation: ln[/(</)] = ln[/(0)] — i? 2 q 2 /3. 
The data reveals that the R s of the SAM-II riboswitch 
changes dramatically with the four sample conditions. In 
the absence of Mg 2+ and the metabolite (condition Al), 
the RNA adopts an extended conformation, which we call 
the unfolded state (U) under low-salt condition (50 mM 
Na + ) with an R g of 3 1 .7 A. In the presence of 10 mM Mg 2+ 
(condition A2), R g decreases to 21.5 A which we designate 
as a compact intermediate state (M). Further compaction 
of R g to 20.7 A occurs in the presence of 0.5 mM metab- 
olite without Mg 2+ (condition A3) to give a conformation 
referred to as the bound state (B). Finally, in the presence 
of both Mg 2+ and the metabolite (condition A4), R g 
compacts even further to 19.5 A which we denote as the 
native folded state (N) (Table 1). Compared with the R g 
value of ~19A calculated from a 2.8-A-resolution X-ray 
crystal structure in metabolite-bound form of SAM-II 
RNA (23), the R g of 31.7 A of the RNA in the absence 
of Mg 2+ and the metabolite is considerably larger, suggest- 
ing that the RNA under this solution condition is 
extended. To further derive the degree of RNA compac- 
tion and highlight RNAs with large unfolded regions, the 
SAXS data were transformed into Kratky representations 
(Figure 4b), where the scattering intensity is weighted by 
the square of the momentum transfer (q) (Ix q 2 versus q). 
According to Porod's law, the scattering profiles from 
folded proteins or nucleic acids with a well-defined 
surface fall approximately as q~ 4 at large q and the 
profiles therefore tend to have a pronounced peak in the 
Kratky representation. In contrast, completely unfolded 
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Figure 4. SAXS analysis of SAM-II RNA under four different states. 

(a) Comparison of experimental scattering profiles of SAM-II. The 
experimental curves were placed on a relative scale in origin. 

(b) Comparison of the SAXS data in Kratky representation for 
SAM-II under the four states, (c) Comparison of P(r) determined 
from experiments using GNOM for the four different states. 



polymers show a linearly rising scattering curve (58). The 
scattering curve of SAM-II RNA under the condition of 
50mM Na-MES buffer without added Mg 2+ or metabolite 
exhibits a flattened plateau with scattering intensity falling 
slightly at large angles. Although it is different from an 
upward rising scattering intensity representative of a com- 
pletely denatured conformation, it is consistent with a par- 
tially 'unfolded state' with an extended conformation. 
This extended conformation is better visualized by 
comparing the distance distribution PDDF P(r), which 
provides information on the shape and maximum intra- 
molecular distance of the particles (58,59). The RNA in U 
state has a remarkably larger maximum intramolecular 
distance f max of 112 A than other states as presented in 
Table 1 and Figure 4c. This D max value, compared with 
the size and the length expected from the SAM-II RNA 
crystal structure further confirms that the RNA conform- 
ation in the U state is not a completely denatured random 
coil, attributed to its partially formed secondary structure. 

Upon addition of 10 mM Mg 2+ to the U state, the RNA 
undergoes significant compaction with a decreased R s 
from 31.7 to 21.5 A and a decreased maximum intramo- 
lecular distance D max from 112 to 72 A, indicating the 
population of a more compact 'intermediate state' 
(M state) that is near folded as confirmed by the appear- 
ance of a concave peak in the Kratky plot. This compact 
M state was also observed in the scattering profile of the 
RNA in the presence of 2mM Mg 2+ without metabolite 
(Supplementary Figures S2 and S3), suggesting that the 
M state is likely fully populated under these solution con- 
ditions and that high stabilizing Mg 2+ ion concentrations 
induce global collapse and structural rearrangement of the 
RNA in the absence of the metabolite. The monovalent 
and divalent ions effectively screen the negatively charged 
backbone of the RNA, facilitating and stabilizing as yet 
unidentified long-range or tertiary contacts within an 
electrostatically relaxed conformation (60-63). 

Remarkably, in the presence of low concentration of 
monovalent ions (50mM Na + ) and 0.5 mM SAM metab- 
olite but in the absence of Mg 2+ ions (B state), the scatte- 
ring pattern is similar to that of the M state, with an 
R g value of 20.7 A and a maximum intramolecular 
distance D mayL of 70 A. However, this B state is distinguish- 
able from those in U and N states, from NMR analysis (see 
below). As indicated earlier, addition of both 0.5 mM me- 
tabolite and 10 mM Mg~ + resulted in further compaction 



Table 1. Overall parameters and quality indicators derived from the scattering data for the SAM-II riboswitch under different conditions 



Structural properties of the RNA under different thermodynamic states 


States 


Solution conditions 


R e (A) 


£>max (A) 


NSD 


Unfolded (U) 


Buffer 


31.7 ± 0.5 


1 12 


0.68 ± 0.01 


Intermediate (M) 


Buffer + 10 mMMg 2+ 


21.5 ± 0.2 


72 


0.62 ± 0.03 


Bound (B) 


Buffer + 0.5 mM SAM 


20.7 ± 0.2 


70 


0.65 ± 0.02 


Native (N) 


Buffer + 10 mM Mg 2+ + 0.5mM SAM 


19.5 ± 0.2 


68 


0.62 ± 0.04 



The R g calculated from the crystal structure (Protein Data Bank ID: 2qwy) is ~ 19 A and the maximum intramolecular distance of the crystal 
structure is 67 A. The R g were fit by Guinier analysis of the low angle scattering data, and the maximum intramolecular distances (£> max ) were 
obtained using the program GNOM. The NSD values of the four low resolution reconstruction models were generated from the program 
DAMAVER. The buffer is 50 mM Na-MES (pH 6.20). 
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of R g to 19.5 A with the smallest D max of 68 A. This value 
of R g is in good agreement with the R g (~19A) calculated 
from the crystal structure (Table 1) (23). Comparison of 
the scattering profiles for the B, M and N states indicates 
that the reduction of R g is accompanied by a small con- 
formational change of the RNA. A sharp peak in the 
Kratky representation of the scattering curve in the N 
state demonstrates that the RNA under this solution con- 
dition is in its native tertiary conformation with a 
well-defined molecular surface (58). 

Low resolution ab initio modeling 

To better visualize the conformational changes, we con- 
structed ab initio models of the SAM-II RNA in the four 
states by using a low resolution structure reconstruction 
method based on the SAXS data using the program 
DAMMIN (see 'Materials and Methods' section) (40). A 
total of 20 independent DAMMIN runs were performed 
for each of the four models corresponding to the four 
states. The resulting 20 models were averaged using the 
program package DAM AVER (41). In this program, the 
NSD between each pair of models was computed. NSD 
provides a quantitative estimate of similarity between 
models: identical models give an NSD of zero, similar 
models give NSD values ml, while models with different 
shapes give NSD values greater than one (64). The final 
low-resolution reconstruction obtained represents the 
average molecular envelope of the ensemble. 

The average NSDs for the four reconstructions are 0.68, 
0.62, 0.65 and 0.62, indicating excellent convergences in 
both individual DAMMIN fits and overall bead model 
ensembles for the RNA under these four different states. 
The low-resolution models corresponding to the four 
states are presented in Figure 5a, which shows different 
views of the final surface of the ab initio models (molecular 
envelope) as mesh representations with 'filtered' models 
rendered as beads and aligned inside the molecular 
envelopes. 

The U state reconstruction reveals an elongated con- 
formation with a -D max of 112 A. Based on the secondary 
structure derived from the X-ray crystal structure, we 
assume that the PI stem is robustly formed in the U 
state RNA, as it is the most established secondary struc- 
ture involving six successive base pairs, and the 3' tail of 
the RNA starting from U31 adopts an open, elongated 
conformation without any secondary or tertiary contacts 
with other parts of the RNA. The PI stem is located near 
the middle section of the reconstruction with a big loop 
from G8 to A24 at one end and the 3' tail at the other. The 
width of the^ middle section of the elongated reconstruc- 
tion is ~20 A, consistent with that of a standard double 
stranded helix from an A form RNA, confirming the as- 
sumption of the existence of a helical stem along the RNA 
molecule. Although it is hard to calculate accurately the 
length of the extended structure, the A-nax of 112 A 
generated from the scattering data is in the right range 
of molecular length. The pairwise NSD values of 
0.66-0.72 of the U state reconstructions indicate that the 
20 models generated from the scattering profile exhibit 
quite similar elongated shape, rather than adopting 



various conformations such as bending or twisting, most 
likely due to electrostatic repulsion under the low salt con- 
dition used here (65). 

The models for the other states exhibit more compact 
conformations with significantly reduced A-nax values 
(72, 70 and 68 A for M, B and N states, respectively), 
and their molecular shapes suggest that Mg 2+ and/or the 
metabolite induced molecular compaction may bring 
about some tertiary interactions or possibly an increased 
extent of secondary structure, as shown by NMR below. 
The N state models show the lowest average NSD value, 
consistent with what we expected to be the most rigid and 
folded surface. Interestingly, the reconstructions from M 
state also show the same lowest average NSD, an indica- 
tion that in the presence of lOmM Mg 2+ the loop region 
and the 3' terminal is highly restrained. The B state recon- 
structions give average NSD value of 0.65, suggestive of a 
larger conformational variation in the absence of Mg 2+ 
than that of M and N states (Figure 5a). 

Comparison between the N state low resolution ab initio 
model and X-ray crystal structure 

To ascertain whether SAXS can provide some tertiary 
structural information, we analyzed the scattering 
profiles from the solution N state and from X-ray struc- 
ture (23). We calculated the scattering data from the 
SAM-II structure determined by X-ray crystallography 
using the program Crysol V. 2.6 and fitted it to the N 
state experimental SAXS curve (Figure 5b) (42). In 
addition, we superimposed the X-ray crystal structure 
(blue ribbon) with the N state ab initio model in surface 
(Figure 5c) representation (the 'filtered' model are 
rendered as beads aligned inside the envelope) using 
SUPCOMB20 (64). The excellent match between the 
curves and structural alignment between the X-ray struc- 
ture and the 'filtered' model indicated that SAXS not only 
provided shape information but also identified the known 
structure in solution. 

We observed a small discrepancy between the model 
and the X-ray structure (Figure 5c). The metabolite 
bound RNA in the crystal structure adopts an even 
more compact conformation than the N state RNA, as 
evidenced by a smaller R g value and a smaller maximum 
intramolecular distance from the X-ray structure 
(Table 1). The tiny divergence might possibly result from 
different experimental conditions. Nonetheless, the RNA 
in solution N state adopts a compact conformation nearly 
identical to that of the X-ray crystal structure. 

Role of Mg 2+ ion in secondary and tertiary structural 
stabilization is pivotal for the RNA to adopt 
well-organized tertiary structure 

Compared with the SAXS data from B state, the slight 
increase of R g and Anax, but otherwise similar scattering 
curve in Kratky plot presented in Figure 4b and global 
shapes shown in Figure 5a, reveals that the conformation 
of the fully populated intermediate M state is nearly iden- 
tical to the bound form, suggesting that either Mg 2+ or the 
metabolite can cause compaction of the overall global 
shape and size of the RNA particle. However, SAXS 
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(a) (b) 




Figure 5. Low-resolution ab initio models of the SAM-II aptamer domain, (a) Average low-resolution ab initio models of the RNA under four 
different folding states: unfolded conformation (first row, blue, U state), conformation in the presence of 10 mM Mg~ + (second row, cyan, M state) 
and in the presence of 0.5 mM SAM metabolite (third row, magenta, B state) and low-resolution reconstruction in the presence of both 10 mM Mg 2+ 
and 0.5 mM SAM metabolite (fourth row, brown, N state). Three different views are shown for the final consensus models (molecular envelope) as 
mesh representations with 'filtered' models rendered as beads aligned inside the molecular envelopes. The black scale bars represent 20 A, the 
diameter of an A-form RNA helix, suggestive of the presence of such helices even in U state, (b) Experimental SAXS profile of SAM-II RNA in 
the presence of the metabolite and lOmM Mg 2+ (cyan) overlaid with the theoretical scattering profile (magenta) calculated from the X-ray crystal 
structure (Protein Data Bank ID: 2QWY). The experimental profile was placed on a relative scale in origin, and the theoretical profile was manually 
adjusted to illustrate the agreement with the experimental profile, (c) X-ray crystal structure (blue ribbon) superimposed with the ab initio model in 
the N state calculated from the experimental SAXS data as surface representation with 'filtered' models rendered as beads aligned inside the 
molecular envelopes. The alignments were performed with SUPCOMB20. 



does not give any detailed indication of the changes in 
secondary or tertiary structure. Therefore, NMR titration 
experiments were performed to determine whether the 
Mg 2+ and/or the metabolite induced structural rearrange- 
ments resulted in any changes of secondary structural 
elements. NMR chemical shift perturbation experiments 
are widely used to map and define binding sites in 
biomolecular complexes, especially in the case of RNA 
for measuring characteristic imino proton and monitoring 
the formation and breakage of Waston-Crick base pairs 
to define conformational change upon interacting with 
metabolites or metal ions (66). ID NMR spectra of the 
titration of the riboswitch with Mg 2+ ion in the absence of 
metabolite reveals that no major changes to the dispersion 
pattern of the imino protons (no chemical shift perturb- 
ation to the existing imino peaks and no strong new 
signals) occur, except of the appearance of one broadened 
but detectable imino signal, even upon addition of 8 mM 
Mg 2+ , indicating that this compaction process does not 
involve any strong base pairing interactions, but results 
from Mg 2+ mediated electrostatic interactions or uniden- 
tified tertiary interactions (Figure 6a). Although the Mg 2+ 



collapsed state has an overall structure similar to those of 
bound form (B state) and native form (N state), it is a 
non-native conformation. 

With titration of the RNA with different stoichiometric 
amounts of the metabolite in the presence of Mg 2+ , 
however, a number of new imino protons peaks charac- 
teristic of canonical or non-canonical Watson-Crick base 
paring appear between 12 and 14.5 ppm, and many pro- 
ton signals resonate around 9 ppm, even upon addition 
of 0.5 mM unsaturating amount of metabolite, 
demonstrating that the metabolite-induced compaction 
in the presence of Mg 2+ indeed involves the formation 
of new base pairing interactions that increases the extent 
of helix packing and that are absent in the M state struc- 
ture (Figure 6b). 

We also carried out NMR titration experiments of the 
RNA with the metabolite in the absence of Mg 2+ , and the 
results again demonstrated that this compaction process 
induced by the metabolite alone is associated with changes 
in base pairing with regards to the appearance of new 
imino signals. Further titration of the RNA with Mg 2+ 
ion in the presence of the metabolite was performed. 
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Figure 6. NMR spectra of SAM-II at 25°C reveal that both Mg 2+ and the metabolite are required to drive folding of the RNA to the native state. 



(a) The imino region of ID H spectra of metabolite-free RNA with different concentrations of Mg 



(b) The imino region of ID H spectra of 



SAM-II riboswithc in the presence of 8mM Mg 2+ with different concentrations of the metabolite, (c) Imino region of the overlaid 2D 'H— 15 N 
HSQC spectra of 13 C, 15 N-uridine-labeled SAM-II riboswitch recorded at 0 or 2mM Mg 2+ concentrations conditions in the presence of 3.2 mM SAM 
ligand. (d) 2D imino 'H- I5 N HSQC spectrum of 13 C, 15 N-fully-labeled SAM-II riboswitch in the presence of 3.2mM SAM ligand and 2mM 
Mg 2+ ion. 



Figure 6c depicts an overlay of two imino 'H— 15 N HSQC 
spectra of 'X/ 15 N-A, U-labeled SAM-II RNA recorded 
at different concentration of Mg 2+ ion. In the absence of 
both Mg 2+ and metabolite, only 1 weak uridine imino 
crosspeak at ~13.4ppm is present (Supplementary 
Figure S4). This is the only uridine present in helical 
stem PI. Note also that all the G-C base pairs in helix 
PI are already formed in the absence of both Mg 2+ ions 
and metabolite. In the absence of Mg 2+ ions, addition of 
the metabolite led to the appearance of five uridine imino 
crosspeaks (blue) at ~11.2, 12.6, 13,1, 13.4 and 14.2ppm 
and caused the disappearance of the weak uridine imino 
peak in the absence of both Mg 2+ and the metabolite, 
suggesting a tremendous conformational change and the 
evolution to a well-folded structure (Figure 6c). Titration 
of the RNA with Mg 2+ in the presence of the metabolite 
shows little change in the guanine imino proton region. 
However, the 'H— 15 N HSQC spectrum in the uridine 
imino region, showed a noticeably change in the 
presence of both Mg 2+ and the metabolite with the 
appearance of one more uridine imino signals at 
~12.8ppm (Figure 6c), indicating that Mg 2+ triggers the 



formation of a new hydrogen bond in the presence of the 
metabolite. 2D NOESY spectrum was recorded for the 
unlabeled RNA in ligand-bound form to perform sequen- 
tial walk (Supplementary Figure S5) and complete the 
assignment of imino protons (Figure 6d). Direct 
hydrogen bonding experiments further confirmed that me- 
tabolite binding triggers formation of new hydrogen 
bonds (Supplementary Figure S6). The imino pro- 
ton sensitive to the presence of Mg 2+ ion is U12 from 
LI loop, demonstrating that Mg~ + is able to cooperatively 
facilitate loop/loop interactions. We suggest that Mg 2+ ion 
binding to U12 likely nucleates the compaction of the 
elongated U state to the compact M state. 

DISCUSSION 

SAXS has proven to be a powerful approach for 
investigating the global size and shape of RNA and, in 
particular, for monitoring the conformational changes in 
the structure of a riboswitch upon binding metabolite 
(67-70). SAXS analysis in combination with NMR spec- 
troscopy enables the resolution of global structures of 
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large RNA and RNA complexes into the underlying sec- 
ondary and tertiary structural constituents (37,59). 
Therefore, we used SAXS to monitor the conformational 
changes of the SAM-II riboswitch in the absence or 
presence of Mg 2+ and the metabolite, and NMR to 
pinpoint elements of structure involved in the global con- 
formational changes. Our SAXS and NMR studies on 
SAM-II riboswitch combined with biochemical results 
suggest a model of RNA folding and functional mechan- 
ism that is characterized by an ensemble of non-native but 
competent conformations for ligand recognition 
(Figure 7). After transcription, the RNA likely adopts 
an elongated conformation (U state, with a R g of 
~32A) at low salt condition (50 mM Na + ). At physio- 
logical condition with Mg 2+ ion in the low millimolar con- 
centration range, the unfolded RNA collapses to the 
M state. 

The RNA in M state samples a restricted ensemble of 
conformations with pairwise NSD values from 0.58 to 
0.71, including conformations resembling the B and N 
states. This inference is supported by the similarity of 
the global shape of the bead model to the two states 
(Figure 5a). The 3' terminus of the RNA folds back to 



form a pseudoknot-like conformation via tertiary contacts 
or Mg + -mediated electrostatic interactions with the LI 
loop region of the RNA. Base pair(s) is (are) formed in 
the process. The exchange between open and transient 
close conformations of the M state RNA, however, is so 
fast that the temporarily formed, dynamical base pair(s) of 
the RNA is broadened as suggested by the induced, weak 
but still visible imino signal at 13.6ppm in Figure 6a. 
SAM-II riboswitch samples different conformations on 
the nanosecond timescale based on a recent molecular 
dynamic simulation study (26). The RNA with this con- 
formation in M state is also somewhat similar to what is 
described as the sensing phase by Stoddard et al. (27), and 
likely monitors the cellular environment for the presence 
of SAM metabolite to adopt a conformation that is favor- 
able for ligand binding. 

A previous smFRET study of SAM-II riboswitch sug- 
gested a stepwise global rearrangement on Mg 2+ ion and 
substrate binding (36). Our SAXS studies now provide 
direct insight into the nature of the global conformational 
changes that accompany these binding events. In 
the presence of low millimolar concentration of Mg 2+ 
ion, the closed conformation becomes increasingly 
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Figure 7. Schematic representation of the regulatory model of the SAM-II riboswitch. Under physiological condition, the transcribed RNA samples 
an ensemble of conformations with its 3' tail folded back to get close to the LI loop via unidentified Mg 2+ -mediated electrostatic interactions. Base 
pairs might be formed between the 3'-tail and LI loop region during the process to form a competent conformation favorable for ligand binding with 
the SD sequence being accessible for ribosome. The open/close process is so dynamic that the formed base pairs are broadened beyond observation. 
In the absence of the ligand, ribosome recognizes and binds the ribosomal binding site. The translation is 'ON'. In the presence of the metabolite, 
base paring of regions (highlighted in green and blue) forms the P2a and P2b secondary structure. Consequently, the SD sequence is sequestered and 
the translation is 'OFF'. 
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populated (36). The M state ab initio model is consistent 
with the state III of the most recent model derived from 
fluorescence studies of SAM-II riboswitch (36) and 
provides direct evidence to support the recently 
proposed mechanism of biomolecular recognition 
designated as 'conformation selection' or 'conformation 
capture' (34,36,71-74). On one hand, in the absence of 
the metabolite, the ribosome can gain access to recognize 
and bind the SD sequence and the translation is turned 
'ON'. On the other hand, in the presence of the metabol- 
ite, ligand binding causes RNA compaction to form a new 
set of base pairs that appear absent or sparsely populated 
in the M state. As a result, the ribosomal binding site is 
occluded and inaccessible to ribosome. The translation is 
turned 'OFF' (Figure 7). The step of ligand binding to 
switch the RNA from M state to N state is believed to 
be rate-limiting since bringing several parts together 
to form a complex binding pocket requires overcoming 
an entropic barrier. Currently, there is no data available 
for determining whether this process is driven thermo- 
dynamically or kinetically. 

Mg 2+ ion induced compaction of the unbound SAM-II 
riboswitch from gel filtration studies on FPLC indicated 
that the Mg 2+ concentration midpoint for the RNA com- 
paction is ~6mM, which is ~2-3x higher than that 
required for forming pseudoknot-like conformation (36). 
In fact, SAXS data was also collected at 2 mM Mg 2+ ion 
concentration, and the ab initio low resolution model is 
quite similar to the M state model, suggesting that the 
unbound SAM-II riboswitch in the presence of 2mM 
Mg 2+ ion also adopts a similar ensemble of conformations 
with pseudoknot-like conformers (Supplementary Figure 
S2). In the presence of the metabolite alone without Mg 2+ 
(B state), the imino 'H- 15 N HSQC spectrum of the RNA 
shows the presence of new signals corresponding to 
uridine imino protons with good spectral dispersion, 
indicating a well-folded conformation. But, compared 
with the spectrum in the presence of both Mg 2+ and the 
metabolite, there are fewer uridine imino peaks indicative 
of less base pairing interactions and suggestive that the 
RNA in B state adopts a relatively less stable conform- 
ation. Even though the RNA has higher degree of protec- 
tion by the newly formed base pairs around the metabolite 
induced binding pocket than that of the unfolded state, it 
retains some flexibility due to the absence of key A-U base 
pairs or Mg 2+ -stabilized electrostatic interactions. 
Surprisingly, the U12 region in LI loop is heavily per- 
turbed by Mg 2+ ion as exhibited by the induced formation 
of new base pairing interaction in the presence of Mg 2+ , 
indicating that Mg 2+ can promote loop/stem tertiary inter- 
actions. This experiment demonstrates unambiguously 
that, in addition to the metabolite, the Mg 2+ ion is also 
required for the RNA to adopt its well-organized native 
folded structure. In the crystal structure of SAM-II 
riboswitch, a magnesium ion is positioned in a pocket 
formed between the 5' termini from two RNA molecules. 
This magnesium ion binding site is far away from the me- 
tabolite binding pocket, suggestive of a delocalized ion 
binding site. However, based on our NMR data, it is 
likely that there is a specific magnesium binding site near 
this triple helix region where the ligand binding pocket is 



formed. Circular dichroism (CD) data also confirmed that 
Mg 2+ is able to improve ligand binding affinity (B. Chen 
and T.K. Dayie, unpublished data). Thus, our studies 
indicate that Mg 2+ and metabolite binding may act 
together to populate a conformation that is competent 
ultimately for turning off translation. 

In the past studies of purine riboswitch aptamer domain 
using smFRET, 2-aminopurine fluorescence, force exten- 
sion experiments, SHAPE, nuclease mapping and NMR 
(28,75-79) arrived at different conclusions about the 
nature of the free and bound states. These differing con- 
clusions likely reflect the use of techniques that probed 
either local or global conformational changes, but not 
both in the same study. Our studies here using both 
SAXS and NMR provide a unified view of the nature of 
the free and bound states of SAM-II that are consistent 
with previous X-ray crystallography, chemical probing 
and smFRET data (23,36). 

In summary, our data demonstrate that the SAM-II 
riboswitch adopts multiple conformations along its 
folding pathway and that under physiological conditions, 
the RNA is unstructured in the absence of the ligand. 
Addition of the metabolite induces the formation of a 
complicated binding pocket with the putative SD 
sequence of the RNA sequestered as well as extensive 
tertiary interactions between the RNA and the metabolite. 
Mg~ + ion is an essential cofactor for proper folding of the 
RNA and can facilitate the formation of a well-organized 
tertiary structure. Occlusion of the SD sequence is 
postulated to inhibit ribosomal binding and to prevent 
translation initiation. Our work complements and 
extends recent studies in two important aspects. First, 
SAXS provides unambiguous evidence for at least four 
conformational states of varying compactions, with 
'Mg 2+ only' induced compaction comparable with both 
'metabolite only' and 'metabolite and Mg 2+ ' induced com- 
pactions; although the underlying secondary and tertiary 
structural arrangements are not completely identical. This 
suggests that the Mg 2+ compacted state samples the 'me- 
tabolite and Mg 2+ induced states in agreement with preva- 
lent ideas about conformational selection as a basis for 
RNA molecular recognition (73,80). NMR analysis 
provides greater insights into the nature of secondary 
and tertiary structural rearrangement of the RNA in the 
presence of Mg 2+ and metabolite, suggesting that the me- 
tabolite binding drives the folded conformer to populate 
the final state competent for shutting down protein 
synthesis. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Figures 1-6. 
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